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C^l ' Gravitational perturbations of neutron stars and black holes are well known sources of gravi- 

» I ' tational radiation. If the compact object is immersed in or endowed with a magnetic field, the 

gravitational perturbations would couple to electromagnetic perturbations and potentially trigger 
synergistic electromagnetic signatures. We present a detailed analytic calculation of the dynamics of 
coupled gravitational and electromagnetic perturbations for both neutron stars and black holes. We 
£S| , discuss the prospects for detecting the electromagnetic waves in these scenarios and the potential 

£N| ■ that these waves have for providing information about their source. 
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I. INTRODUCTION 



Multi-messenger astronomy has arrived. Already astro-particle observations (neutrinos and cosmic rays) are comple- 
menting traditional electromagnetic observations. The third pillar is almost ready with near future gravitational- wave 
observations by interferometric detectors like LIGO, Virgo, GEO600 and LCGT [H, 0. This new astronomy will en- 
able multi-channel observations of astrophysical phenomena such as 7-ray bursts, supernovae, or flaring magnetars, 
unveiling an unprecedented view of the nature of the source and its environment. 

An important component in many astrophysical phenomena is strong magnetic fields, as demonstrated by the active 
role they play in the accretion processes of low-mass X-ray binaries and GRBs Q . The presence of strong magnetic 
fields opens up the possibility for interesting effects. Among them, which is the central topic of this work, is the coupling 
' between electromagnetic and gravitational emissions that could yield synergistic multi-messenger observations. In 
, particular, it is important to assess the conditions in which electromagnetic and gravitational emissions influence 
■ each other. There are already hints for such scenario. It is believed that the flare activity of magnetars seems to be 
(T) \ associated with starquakes [4| ■ These quakes are responsible not only for dramatic perturbations and rearrangements 
1 of the magnetic field, but also for the breaking of the neutron star crust and internal motions, possibly resulting in 
the emission of gravitational waves. Detailed studies of magnetar flare activity have revealed a number of features in 
the afterglow, which can be associated with the crust oscillations as well as with Alfvcn waves propagating from the 
core towards the surface [5l-fl7l|. 

The link or coupling between electromagnetic radiation and gravitational waves have been investigated for some 
cases. One of them looked at the propagation of gravitational waves linearly coupled to an external magnetic field fl8j . 
It was shown that this configuration triggers magneto-hydrodynamics waves in the plasma fl9l - l23j . Furthermore, the 
linear nature of the coupling limits the electromagnetic waves to low frequencies, in the best case a few tenths of 
kHz, which will be easily absorbed by the interstellar medium or plasma. In order to produce high frequency and 
detectable electromagnetic waves, non-linear couplings are needed, requiring much stronger gravitational waves. In 
most of these studies, the gravitational waves were assumed to propagate on a flat space-time background. This is 
a reasonable assumption when the interaction between the gravitational and electromagnetic waves takes place far 
from the source. There have been only very few attempts to treat the electromagnetic-gravity coupling in the strong 
field regime [Hj|. 

The aim of this work is to study the interaction of electromagnetic and gravitational waves in the vicinity of mag- 
netized neutron stars or black holes immersed in strong magnetic fields using perturbation theory, paying particular 
attention to how gravitational modes drive the excitation of electromagnetic perturbations. Our work also includes 
estimates of the energy transferred between the gravitational and electromagnetic sectors. As expected, we find that 
the excited electromagnetic waves have roughly the same frequency as the driving gravitational waves, i.e., of the 
order of a few kHz. Electromagnetic waves at these low frequencies can be easily absorbed by the interstellar medium. 
As a consequence, one needs to associate them with secondary emission mechanisms (e.g., synchrotron radiation) in 
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order to be able to trace the effects of gravitational waves on the strong magnetic fields. The later process can be 
studied following the mechanisms described in [l9l - l23j . and there is work in progress for the special case of strong 
gravitational fields. 

This article is organized as follows: Section [TT1 gives details of the space-time background configuration. In Sec. IIII1 
we review the general form of the perturbation equations, their couplings, and the angular dependences of the various 
types of electromagnetic and gravitational perturbations. In Sec. IIV| we reduce the equations to the particular case 
of dipolc electromagnetic perturbations driven by the quadrupole gravitational mode for the case of a neutron star 
background. In Sec. [V] we do the same as in Sec. IIVI but for the case of a black hole and consider both the case of 
axial and polar gravitational perturbations. In Sec. IVII we show numerical results or dipole electromagnetic waves 
driven by quadrupole gravitational waves with axial parity for both neutron stars and black holes. Conclusions are 
given in Sec. lVIII We adopt geometric units, c = G = 1, where c and G denote the speed of light and the gravitational 
constant, respectively, and the metric signature is (—,+,+,+). 



II. EQUATIONS FOR THE BACKGROUND 

The background space-times we are considering (neutron stars and black holes) are governed by the Einstein- 
Maxwell equations, which read: 

Guv = Sn (T^ + E^) , (2.1) 

(T^ + E^). v = 0, (2.2) 

F» v . v = 4ttJ^, (2.3) 

F^v,\ + F\n,v + Fv\,n = , (2.4) 

The tensors that appear in these equations are: The Einstein tensor G^, the Faraday antisymmetric tensor F llVl the 
electromagnetic four-current J 1 *, the energy- momentum tensor of the matter fluid and the energy-momentum 
tensor of the electromagnetic field is E^ v . The energy- momentum tensors are explicitly given by 

T MV = (/£> + p)u^u v + pg^ , (2.5) 

E^u = ^ [g pa F w F av - -^F^FP^j , (2.6) 

where p stands for the energy-density, p for the pressure, and it M for the four-velocity of the matter fluid. 

The presence of a magnetic field could in principle induced deformations to the neutron star or black hole we are 
considering. However, even for astrophysically strong magnetic fields, B ~ 10 16 G, as in the case of magnctars, the 
energy of the magnetic field £b is much smaller than the gravitational energy £q, by several orders of magnitude. 
In fact, £b/£g ~ 10~ 4 (£>/10 16 [G]) 2 . Therefore, in setting up the background space-time metric, one can ignore the 
magnetic field. That is, the background metric has the form 

ds 2 = -e v di 2 + e x dr 2 + r 2 (d0 2 + sin 2 9d(/> 2 ) , (2.7) 

where the functions v(r) and A(r), in the interior of a neutron star, are determined by the well-known Tolman- 
Oppenheimer-Volkoff (TOV) equations (see, e.g. (25|) and the matter fluid four-velocity = {e~ v / 2 , 0, 0, 0). In the 
exterior of a neutron star, and in the case of a black hole, they are determined by the standard Schwarzschild solution: 
er x = e u = 1 - 2M/r. 



A. A Dipole Background Magnetic Field: Exterior region 



Next, we compute the magnetic field for both the neutron star and the black hole. We consider first the exterior 
(vacuum) solution. In this case, the component of Maxwell equations given by Eq. (|2.4p is automatically satisfied. 
The magnetic field is then obtained by solving the remaining Maxwell equations, Eqs. (|2.3[) . which in vacuum reads 

F> iV . v = Q, (2.8) 

with F^ = A VtlJt — A^ >v . Since the background space-time is static, it is natural to assume that the magnetic field is 
also static. In addition, we require the magnetic field to be axisymmetric and poloidal, 

= (O, e- A /2 B (-) ( r ) cos 9i e -V2 s M ( r ) sin 9j ) , (2.9) 
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which has a dependence on the polar coordinate, 6. From the relation between the magnetic field, the matter fluid 
velocity u M , and the field strength 



(2.10) 



where t^vaf} is the complete antisymmetric tensor, determined by the convention eoi23 = ^ ^ s n0 ^ difficult to 

show that the only non- vanishing component of the vector potential is the 0-component, which we will denote as 
^(ox) therefore, the vacuum Maxwell equation (|2.8|) in the Schwarzschild background becomes 
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2M\ dA 
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Expanding A^?* 1 in vector spherical harmonics as 



A™ =a\™\r)Bm0d Pl M (coB0), 



we rewrite Eq. (|2.1ip as 



dr 



2M \ da] 
r J dr 



(ex) 



f(Z+l)a£>=0. 



The solution of this equation for the dipole case (Im = 1) has the form [26 
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where fid is the magnetic dipole moment for an observer at infinity. With the solution of Eq. (|2.14[) and Eq. (|2.10l) , 
the coefficients of the magnetic field in Eq. (|2.9p are given by: 



B[ ox \r) 



Bt\r) 



Notice that in the limit r — > oo, 



2a 



(ex) 



3Md 


In 


\-™) 

\ r J 


2M 
H 


2M 2 ~ 


4Af 3 


r 




3/Ud 


■„| 




M 


M 


4M 3 r 


H + 

r 


r - 2M 



B[ ex \r) 



2fi a 



and B ( ™\r) 



,•4 



(2.15) 
(2.16) 

(2.17) 



B. A Dipole Background Magnetic Field: Interior region 

We assume that the magnetic field inside the star is also axisymmetric and poloidal, with current J ;t = 
(0, 0,0, Jcj,) [13, H3|- The ideal MHD approximation is also adopted, i.e. infinite conductivity er, which leads to 
Efi = F^u 1 " = 0, as follows from the relativistic Ohm's law 

4-7T 

Ffiuu" = — ( J M + u„J u u v ) . (2.18) 
a 

Therefore, the vector potential A^ is similar to that for the exterior magnetic field, i.e. A^ — (0, 0, 0, )• The 
counterpart equation to Eq. (|2 . 1 1 1) but for the interior is 



e 



-X 



8 2 A^ ] ld 2 A { r } e~ x 8A^ lcos05A (in) 



+ V _ A ~H 51 3—77^ = ~**J* ■ (2.19) 



dr 2 r 2 d6 2 2 dr r 2 s'md 86 



Expanding both, the vector potential A^' and the current J^, in vector spherical harmonics, one gets 



A^\r,6) = 4™\r) S m6deP lM (co S e), (2.20) 
Mr,0) = j iM (r) sm6 d e P lM (cos6), (2.21) 
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which can be use to rewrite Eq. (|2.19[) as 



,(in) 



dr 2 



dr 
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(2.22) 



It is only feasible to obtain numerical solutions to Eq. (|2.22p . even for the dipole case (Im = 1), since among other 
things the coefficients are also computed numerically from the TOV equations. In addition, when prescribin g j\ (r), 
it must satisfy an integrability condition (see [2^, [30[ for details). We adopt a current with a functional form [281] : 



3l(r) = for 2 (p + p) , 



(2.23) 



where /o is an arbitrary constant. In addition, we should impose the following regularity condition at center of the 
neutron star, 



<4 in) =a c r 2 +0(r 4 ), 



(2.24) 



where a c is also an arbitrary constant. These arbitrary constants, /o and a c , are determined by from the matching 
conditions at the surface of the star, namely that a± and ai, r are continuous across the stellar surface. Finally, once 
we have the numerical solution for a\ (r) , the magnetic field is obtained from 



= (o,e- A / 2 Bf n) (r)cos0,e- A / 2 ^ in) (r)sin0,o) 



with 



9 ( in ) 

B[ in \r) = 



and B { 2 in) (r) = - 



(2.25) 



(2.26) 



With the magnetic field determined both in the interior and exterior regions, the Faraday tensor for the background 
field becomes 



F„ 



B a i 



r 2 sin ( 



/ 



\0 -B 2 sin6» Bicosfl 




(2.27) 



III. PERTURBATION EQUATIONS 



We consider small perturbations of both the gravitational and electromagnetic fields, which can be described as 



F, 



= Fau + / 



(3-1) 
(3.2) 



where g^ and F^ v are the background quantities derived in the previous section. The tensors and denote 
small perturbations, i.e. h^ u — 8g^ v and f^ v = SF^ IU . Linearization of the Einstein-Maxwell equations yields 



SG 



S (T M f v + E^ 



(-g) ll2 r 



8ttS (T^+Epv) , 

o, 

AttS \(-g) 1/2 jH - d v \F» v 6(-g)V 2 
0. 



(3.3) 
(3.4) 

(3.5) 
(3.6) 



From Eq. (|3.5I) . we find that the electromagnetic perturbations are driven by the gravitational perturbations via 
the term containing S(-g) 1 ^ 2 in the right hand side. On the other hand, for simplicity, we omit the back reaction of 
the electromagnetic perturbations on the gravitational perturbations, i.e. we set SE^ = 8{E^ V . V ) = in Eqs. (|3.3j) 
and (|3.4p . This simplification is based on the assumption that the energy stored in gravitational perturbations is 
considerably larger than that in electromagnetic perturbations, which are typically driven by the former. On the 
other hand, in the giant flares of SGR 1806-20 and SGR 1900+14 [3lT[33l ]. whose peak luminosities are in the range of 
10 44 — 10 46 ergs s _1 , the dramatic rearrangement of the magnetic field might lead to emission of gravitational waves. 



Nevertheless, recent non-linear MHD simulations [34 38] do not support these expectations. 
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The first two perturbativc equations, Eq. (|3.3[) and Eq. (|3.4[) . have been studied extensively in the past, in the 
absence of magnetic fields, both for stellar and black hole backgrounds (see, e.g. [3^-113 )• Thus, in this article we use 
the perturbation equations derived in earlier works, and we derive the analytic form of the perturbation equations for 
the electromagnetic field together with their coupling to the grav itational perturbations. 

The metric perturbations , in the Regge- Wheeler gauge [39( , can be decomposed into tensor spherical harmonics 
in the following way 
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(3.7) 



where -ffo.im, H\ t i mi H2,i m and K\ m are the functions of (i, r) describing the polar perturbations, while /io,im and fti,im 
describe the axial ones. On the other hand, the tensor harmonic expansion of the electromagnetic perturbations, f^u, 
for the Magnetic multipoles (or axial parity) are given by 
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while the expansion for the Electric multipoles (or polar parity) can be written as 
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(3.9) 



Hereafter, the quantities describing the magnetic- and electric-type electromagnetic perturbations will be denoted 

with the indices (M) and (E), respectively. We point out that h^ v is a symmetric tensor, while both ffffi and ff$ 

are anti-symmetric tensors, i.e. ff£P = —fvft and fjb/ = —fiffi. From the perturbed Maxwell equations, Eqs. f|3.6j) . 
we can obtain the following relations connecting the above perturbative functions: 
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(3.10) 
(3.11) 



where A = 1(1 + 1) 



Notice that f^} m 



and VP, defined as 



Lf( E ) 
^ ■'01, Im ' 



(3.12) 



are gauge invariant variables (see Eq. (11-27) in Ref. [45[ and Eq. (11-11) in Ref. [4f 



A. Perturbations of a Dipole Magnetic Field: Exterior region 



In the exterior vacuum region, we adopt the condition S J M = 0. With this condition, the perturbed electromagnetic 
fields will be determined via the linearized form of Maxwell's equations, Eqs. p.5[) . (assuming that J' 1 = 5J^ = 0) 



(-g) l,2 r 



{-g) 1/2 F^g^h afi 



(3.13) 



together with the perturbed Maxwell equation Q3.6P . Equation (|3.13[) for fi = t and fj, = r can be written down as 

E [A^Yi™ + A\ m A) cos9Y lm + smed e Y lm + C^d^} =0 (7 = 0, 1) , (3.14) 



l,m 
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where the indices "^4" and "P" stand for axial and polar gravitational perturbative quantities, and obviously "7" 
stands for the t and r components of Eq. p,13[) . The coefficients of Eq. (|3.14p have the following expressions 
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One can decompose the equations above for a specific mode with fixed harmonic numbers (i,m), by multiplying with 
Yj^ and integrating over the two-sphere, i.e. 
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In a similar way, from the two remaining equations, i.e. Eq. f|3 . 13[) for (i — 9 and fi — <f>, one gets the relations 
^ | (aim + ai m cos 9) deYim - U3i m + /3; m cos 6>) (d^Yim/ sin 0) + rji m sin 0Y" ;m + xzm sin 9Wi m j = 0, 
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These equations lead to an extra set of evolution equations for a specific mode (I, m) by multiplying with Y* m and 
integrating over the two-sphere: 
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where the coefficients are given by 



a, --(\'-v')f {E) + e A " l/ f (E) - f (E) ' 
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B. Perturbations of a Dipole Magnetic Field: Interior region 

In the stellar interior, because we have adopted the ideal MHD approximation for which F^vu" = 0, the components 
of the perturbed electromagnetic field tensor are determined by using the perturbed Maxwell equation (|3.6[) . i.e. 



fo„ = e»l 2 F, LU 5u» , 
where 6u^ is the perturbed fluid 4-velocity, defined as 

Su* = ^e~ u/2 H ,i m , Rim, Vimde - U im sin" 1 6»<9 , V hn sin" 2 (9<9 + U hn sin" 1 6d e 

From Eq. (|3.37[) one can get the following equations 



Yl {fo?i,n Y im - r 2 B 2 e v ' 2 {V tm d^Y lm + U lm Bm0d Y lm )} = , 
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(3.37) 
(3.38) 

(3.39) 
(3.40) 
(3.41) 



where the coefficients Aim and Bi m are functions of the perturbed electromagnetic fields, while Aim, Blm, and Ch 
are functions of the perturbed matter fluid 4-velocity. The expressions for these coefficients are 



f (E) 
J02,lm ' 



AM) 

hn.im i 



Aim 

Blm = 

Aim = r 2 B ie ^ 2 Ulm 

Blm = -^B^^Vlr 

Cim = r 2 B 2 e v/2 Rim 



(3.42) 
(3.43) 
(3.44) 
(3.45) 
(3.46) 



By multiplying Eqs. (|3~39)) . (|3~40| . and ([3~4l"j) with Y* m and integrating over the two-sphere we can obtain the following 
system of equations that depends only on r 

/oSL ~ r 2 B 2 e^ 2 [imVim + Q tm (l - l)Ui- lm - Qi+im(l + 2)U l+lm ] = , (3.47) 
KA im[Bi m + Ci 

tn\ Qlm 

(l - 1)(Z + l)Al- lm + Ql+lmlil + 2)Ai +lm = , (3.48) 

ABlm + imAlm + Qlm{l + ~ l)Bl-lm ~ Cj-lm] + Ql + lrJ[(l + 2)Bl + lm + Cl + lm] = , (3.49) 
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where 

/ (I - m)(l + m) 
(2Z-1)(2Z + 1) 



lira 



(3.50) 



Finally, we should compute Eqs. (pT23j) . (pT28|) . and (pT29j) for the exterior region, and Eqs. ([3^47]) . (f3T48|) . and (pT49j) 
for the interior region of the star. From this system of equations, we can see the specific couplings between the 
electromagnetic and gravitational perturbations. For example, an electromagnetic perturbation of specific parity 
with harmonic indices (Z, m) depends on the gravitational perturbations of the same parity with (I, m) as well as the 
gravitational perturbations of the opposite parity with (Z ± 1, m). In other words, for the special and simpler case of 
axisymmetric perturbations (m = 0), we arrive at the following conclusions: 1) Dipole electric (polar) electromag- 
netic perturbations will be driven by axial quadrupole gravitational perturbations, and 2) Dipole magnetic (axial) 
electromagnetic perturbations will be driven by polar quadrupole and radial gravitational perturbations. These two 
types of couplings will be discussed in detail in the next sections. 



C. Junction conditions for perturbed electro-magnetic fields 

In order to close the system of equations derived in the previous subsection, we should impose appropriate junction 
conditions on the stellar surface. Such junction conditions for the perturbed electromagnetic fields can be derived 
from the conditions 

„/*£B(i») = n^SB^, (3.51) 
qfSEW = q M v 6E<r> , (3.52) 

where n M is the unit outward normal vector to the stellar surface, while q^ is the corresponding projection tensor 
associated with . These junction conditions lead to the following set of equations: 

j(M)(in) _ ^(M)(ex) ^ ^ 

f MM = /(f)C-)= 0, (3.54) 
/£P (ex) = 0. (3.55) 



IV. DIPOLE PERTURBATIONS OF A MAGNETIC FIELD ON A STELLAR BACKGROUND 

In the previous section, we provided the general form of the perturbative equations. In order to focus on a simple 
case, we only consider axisymmetric perturbations (m — 0) in this section. In this way, the various couplings become 
less complicated. Under these conditions, we study the excitation of dipole electric perturbations driven by axial 
gravitational ones and dipole magnetic perturbations driven by polar gravitational ones. These perturbative modes 
are actually the most important ones from the energetic point of view. 



A. Dipole Electric Perturbations driven by Axial Gravitational Perturbations 

Here, we consider only dipole "electric type" perturbations driven by quadrupole axial gravitational perturbations. 
Since we neglect the back reaction of electromagnetic perturbations on the gravitational ones, the quadrupole axial 
gravitational perturbations of a spherically symmetric star can be described by a single wave equation [4fJ, |47| , which 
is given by 

d 2 X lm d 2 X lm v (k 6m A 

where 

X lm = h hlm and il= e ^)/ 2 ^. (4.2) 

r or or* 
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Note that r* is the tortoise coordinate defined as r* = r + 2M \n(r /2M — 1). Since there are no fluid oscillations if the 
matter is assumed to be described as a perfect fluid (unless we introduce rotation), the spacetime only contains pure 
spacetime modes, i.e. the so-called w- modes [IH, E3, EH . In this case, the axial component of the fluid perturbation, 
Ui m , has the form Ui m = —e~ 1/ ' 2 hoj m /r 2 , while the component of /io,im is computed from the equation 



d_ 

Of 



a 



"Mm, + r- Xi m , 



(4.3) 



which is used later to simplify the coupling terms between the two types of perturbations. 

On the other hand, in the same way as in the case of electromagnetic perturbations in the exterior region, Eqs. 

()3.1ip and (|3.23|) for I = 1, and Eq. (|3.28j) lead to three simple evolution equations for the three perturbation functions 

f (E) f (E) and f (E) • 
./ 12,10' J01,10' dIiu •> 02,10" 



a /l2,l( 

at 



02,10 



dr* 



f (E) 
^01,10 ' 



10 



at 



re"/: 



(E) 
12,01 



s. 



(1) 

20 ' 



(E) 
02,10 



at 



, U J 12,10 

<9r* 



4- i/p 2 » f (E) 
+ v e Jl2,10 



,(2) 



(4.4) 
(4.5) 
(4.6) 



where S 20 ^ and S%q are the source terms describing the coupling of the electromagnetic perturbations with the 
gravitational ones, and are given by 



c(i) 

°20 



3Q 2 



-Bi - e"B 2 X 



^20 



rB- 



dX 2 



S. 



2) 



-Q2ore" Bi X 20 



(4.7) 
(4.8) 



In order to derive second-order wave-type equations for the electromagnetic perturbations, we introduce a new func- 
tion: * im = *j m (i,r), given by 



L J 121m 



With this variable, the above evolution equations can be written as 



9* 



10 



at 

u J 01.10 



df (E) 

Lf iQ2,10 

dr* 



e JOl.lm ' 



at 

■(E) 
02,10 



r*1 



S. 



a/ ( (E) 



r 



10 



5, 



20 ' 



(2) 



(4.9) 

(4.10) 
(4.11) 
(4.12) 



dt dr* 

From this system of evolution equations, one can construct a single wave-type equation for the "electric" perturbations 



d 2 ^io 
dt 2 



drl 



-n e Wio — 5 20 , 



(4.13) 



where the source term S 20 is given by 



? (E) 
'20 



as. 



(2) 



dr* 



e v S, 



(D 

20 



(4.14) 



Without the coupling term, this wave equation outside the star is the well-known Rcggc- Wheeler equation for elec- 
tromagnetic perturbations. It should be pointed out that ^ is not a gauge-invariant quantity while the function \f 
given by Eq. (|3 . 1 2[) is a gauge invariant variable, where both variables \& and ^> can be related to each other via the 
evolution equation (|4.1ip . i.e. 



(4.15) 



10 



Finally, the electromagnetic perturbations in the interior region are determined from Eqs. (|3.47l) . (|3.48[) . and (|3 . 1 1 1) . 



i.e. 



f( E ) _ R c(3) 
J01,10 — D ? D 20 ' 

f( E ) - Ir cj(3) 

^02,10 ~~ 2 1 20 ' 
L, J12,10 _ "J 02,10 j.(E) 

Joi,io i 



9t 



<9r 



where 



4? = -3Q 20 r 2 e^ 2 f/ 2 



(4.16) 
(4.17) 

(4.18) 
(4.19) 



B. Dipole Magnetic Perturbations driven by Polar Gravitational Perturbations 



As it was mentioned earlier in Sec. IIIII for the case of axisymmetric perturbations, the "magnetic (axial) type" 
perturbations of the electromagnetic field with harmonic index I are driven by polar gravitational perturbations with 
harmonic index Z±l. Here, we consider the axisymmetric perturbations (to = 0) for the dipole (I = 1) electromagnetic 
fields, which are driven by quadrupole (I = 2) gravitational perturbations. 

For the descri ptio n of the perturbations of the spacetime and the stellar fluid, we adopt the formalism derived by 



Allen et al. in [44|. In this formalism, the perturbations are described by three coupled wave-type equations, in 



such a way that two equations describe the perturbations of the spacetime and the other one the fluid perturbations. 
In addition to these three wave equations, there is also a constraint equation. The two wave-type equations for the 
spacetime variables are 



dt 2 



dr 2 



2e v 4e 2l/ 
— [2Trr 3 {p + 3p) +m-(n + l)r] S hn = — — 



(m + Airpr 
r — 2m 



3\2 



4irpr — 3m 



Ft 



lm 



(4.20) 



lm 



d 2 F lm 2e v 



dt 2 



dr 2 



— [27rr 3 (3p + p) + m - (n + l)r] F t 



I in 



2 [4nr 2 (p + p) - e~ x ] S lm + 8n(p + p)re» (\ - 



where F im , Si m , and Hi m are given by 



Fi m (t,r) = rK Wl 



Si m (t,r) = — (H A m ~ Ki m ) 

r 



Hi m (t,r) 



5p 



lm 



P + P 



(4.21) 

(4.22) 
(4.23) 

(4.24) 



while 6pim is the perturbation in the pressure, n = (I — l)(l + 2)/2, and C s is the sound speed. On the other hand, 
the wave equation for the perturbed rclativistic enthalpy Hi m , describing the fluid perturbations, is 



1 d 2 H lr 
C? dt 2 



d 2 H ln 
dr 2 



,(i/+A)/2 



(m + Airpr 3 ) ( 1 



1 

C 2 



+ 2(r - 2m) 



dHu 

dr* 



2e h 



2nr 2 (p + p)[:i+—)-{n + l) 



(m + Airpr A ) ( 1 - — 



1 \ e^- v )l 2 ( e v OF, 



2r 



r 2 dr* 



dr* 



(to + Airpr 3 ) 2 
r 2 (r — 2m) 



Tpr3 f i-&) 



2r 2 



2(to + Airpr 3 ) 2 1 to -+- Atrpr 3 ( 1 

1 ~ 7^? I ~ 47r7 W + P) 



r 2 (r — 2m) C 2 



2r 2 



C 2 



F,, 



(4.25) 
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This third wave equation (|4.25j) is valid only inside the star, while the first two are simplified considerably outside 
the star, which can be reduced to a single wave-type equation, i.e. the Zerilli equation (see 44| and §VB|. Finally, 
the Hamiltonian constraint, 



d 2 F ln 
dr 2 



—2 (m + 4irr 3 p) — — + -3 [l27rr 3 p - m - 2(n + l)r] F lr 



dr„ 



4m 



8Trr\p + p) - (n + 3) H 



87T7* 

Sim + ^e v (p + p)H lrn = , (4.26) 



can be used for setting up initial data and monitoring the evolution of the coupled system. 

Regarding the quadrupole gravitational perturbations, the perturbation equation for the "magnetic type" dipole in 
the exterior region is obtained from Eq. p.29[) as 



d 2 <S> lt 

dt 2 



d 2 <S> 



10 



dr 2 



10 



? (M) 
'20 ' 



(4.27) 



where <f> 



,(M) 

'23,im 



and 



S. 



(M) 
20 



(2B 2 + rB 2 ') r 2 S 2 



e v B 2 + re v B 2 ' - -B 1 F- 



20 



rB 2 



dF 



■211 



dr* 



(4.28) 



In order to derive the wave equation (|4.27[) . we have used Eq. (|3.10p and the (r, </>)-component of the perturbed Einstein 
equations, i.e. e~ v Fl\ — Hq + K' — v' Hn = 0. We remark that the wave equation (|4.27[) without the source terms 
is the same as the one derived in [4j| In addition, the other components of the electromagnetic perturbations, 
/1210 an d /021O' can ^ e determined with $10 via the relation (|3.10[) . 

Finally, from Eq. (|3 .49[) and Eq. (|3 . 10[) . we can obtain the equation that determines the dipole "magnetic type" 
perturbations for the interior region: 



<9$ 



10 



dt 



Q 20 r 2 e»/ 2 (B 2 R 20 - WMo) , 



where the perturbations of the fluid velocity, R 2 q and V20, in the source term are given by 



OR 



dt 

dV 20 
dt 



e u/2-\ 



e~ v S 2 



2r 2e 



>'/- 



lip + 3p Zrv' 
2(p + p) + ~ 

re 620 H 2H 20 

p + p r 



3p- P 
2r 2 (p + p) 



(F20 - rF^o 



H 20 



(4.29) 

(4.30) 
(4.31) 



V. PERTURBATIONS OF DIPOLE MAGNETIC FIELD ON A BH BACKGROUND 

The perturbations of a dipole magnetic field on a Schwarzschild black hole background are described by the same 
set of perturbation equations as in the exterior region of the star except for the boundary conditions, i.e. the 
boundary conditions for the neutron star imposed on the stellar surface are Eqs. (|3.53[) - (|3.55p . while for the 
black hole case one should impose the pure ingoing wave conditions at the event horizon. Then, even in the case 
of the black hole background, we observe the same coupling of the various harmonics of the electromagnetic and 
gravitational perturbations as for the neutron star background. That is, for the axisymmetric perturbations, the 
"electric" dipole (1 = 1) perturbations of the electromagnetic fields will be driven by axial quadrupole (/ = 2) 
gravitational perturbations, while the "magnetic" dipole (I = 1) perturbations of the electromagnetic fields will be 
driven by polar quadrupole (I = 2) gravitational ones. In this specific case, our study is similar to the work in [24j . 
although they use a different formalism. 

A. Dipole Electric Perturbations driven by Axial Gravitational Perturbations (BH) 

The axial quadrupole (/ = 2) gravitational perturbations are described by the Rcggc- Wheeler equation 

X lm = 0, (5.1) 



d 2 X lm d 2 X lm , „ f A 6M 



dt 2 dr 2 
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where 

e v 

Xim = — hum . (5-2) 

r 

In accordance with the results of Section IIV Al the perturbations of the electromagnetic fields will be described by a 
single wave equation, that is, the Regge- Wheeler equation for electromagnetic perturbations, give by 

1 -^e *i = SV , (5-3) 



dt 2 dr 2 r 2 



where the source term becomes of the same form as in Section [TV At ^>i m = e v 



lm ' 



B. Dipole Magnetic Perturbations driven by Polar Gravitational Perturbations (BH) 

The equation describing the "magnetic" type perturbations driven by the gravitational perturbations is the same 
equation as the one derived for a neutron star background (see Eq. (|4.27p ~). that is 



~dt 2 drT + ^ e ' <Sw - b ^ 



3e'*io = S& V , (5.4) 



where $/ m = f%fi m , and the source term is also of the same form as in Eq. (|4.28|) . The perturbative equation for the 
spacetime variables can be written in the form of the Zerilli equation 

^-^ + Vz(r)Z lm =0, (5.5) 

. , 2e» [A 2 (Ai + l)r 3 + MIA{r 2 + 9Al 2 A ir + 9M 3 1 
Vz{r) = r 3 (rA 1+ 3Af) 2 ' (5 - 6) 

where Ai = (I + 2) (I — l)/2. Meanwhile, in the same way as for the neutron star background, one can also adopt Fi m 
and Si m as the perturbation variables for the spacetime. In this case, the two wave equations simplify to become 



d 2 Sl m 


9 2 Slm „ | 

~ drl +e 1 


(A 


dt 2 


\r 2 


d 2 Fim 


d 2 F lm + e „ 


fA 


dt 2 


dr 2 


^r 2 



A 2M 



Sim = -— e v I 3 - — Fim , (5.7) 
Fim = -2e v S lm , (5.8) 



which have to be supplemented with the Hamiltonian constraint equation 

d 2 F lm MdF lm A dS lm 1,,„, AW n 

—^-n 5 7{ o e F lm-r— (4e + A) S hn = . (5.9) 

c/rj r z or* r z or, 2 

Note that there are useful relations between the perturbation variables {Fi m , Si m ) and the Zerilli function (Z), i.e. 

dZim , Ai(Ai + l)r 2 + 3AiMr + 6M 2 

Fim = r ^7 + r(A ir + m Zlm ' (5 - 10) 

IdFim {Ai + 2)r - M (A 1 + l)(A 1 r + 3Af) y 

*im — ; 3 firnl 3 ^Zrn , l- 3 - 11 ] 

which can be used in constructing initial data (since the Zerilli function is gauge invariant and unconstrained) , or for 
the extraction of the Zerilli function. 



VI. APPLICATIONS 



As an application, we consider the case in which dipole "electric type" perturbations are driven by axial gravitational 
ones and present numerical results. First, we study the coupling on a Schwarzschild black hole background and later 
on the background of spherical neutron stars, as discussed in W Al and in £ )IV Al respectively. The more complicate 
cases that involve the driving of "magnetic type" electromagnetic field perturbations driven by polar gravitational 
ones will be discussed elsewhere in the future. 
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A. Perturbations on a Black-Hole Background 



In order to calculate the waveforms in the black hole background, we need to modify the background magnetic field 
near the event horizon. The reason for this is that the solution for a dipolc magnetic field in vacuum diverges at the 
event horizon (see Eqs. (|2.15l) and (|2.16[) h In fact, the isolated black hole cannot have magnetic fields due to the no 
hair theorem. But, according to the simulations of the accretion onto the black hole, the magnetic field can reach 
almost to the event horizon, because the accreting matter will fall into the black hole with infinite time [5l|, [52|. Thus, 
we adopt a simple modification of the dipole magnetic field near the event horizon, that is, we set B±(r) = Bi(6M) 
and -62(f) = I?2(6M) for r < 6M, where the position at r = 6M corresponds to the innermost stable circular orbit for 
a test particle around the Schwarzschild black hole. The magnetic dipole moment fid is identified with the normalized 
magnetic field strength B15, defined as S15 = B p / (10 15 [G]), where B p is the field strength at r = 6M and 6 = 0. 

We assume vanishing electromagnetic perturbations i.e., \I , io = d^io/dt = at the initial time slice t = 0, 
while the initial gravitational perturbations, X20, are prescribed in terms of a Gaussian wave packet. Under these 
initial conditions, the electromagnetic waves will result from the coupling to the gravitational ones. In the numerical 
calculations, we adopt the iterated Crank-Nicholson method [53[ with a grid choice of Ar* = Q.1M and At = Ar»/2 
(see [54| for the dependence of the choice of Ar* and At on the waveforms). 

The energy emitted in the form of either gravitational (Egw) ° r electromagnetic waves (Eem), is estimated by 
integrating the luminosity ; ) for the axial gravitational waves and for electric type electromagnetic waves (i^i t ) > 



which are described by the following formulae [45l |46| 

(A) 1 (J -2)! 



L 
L 



GW ' 1 16tt(Z + 2)! 
e) 1 (1 + 1)! 



dX, 



10 



EM,! 



4tt (I - 1)! 



dt 



dt 



(6.1) 
(6.2) 



In practice, we can find a relation between the energy emitted in gravitational and electromagnetic waves for a given 
initial spacetime perturbation, which has the form: 

Eem = azBi 5 2 Ecw , (6.3) 

where a is a "proportionality constant" . 

In the simulation that we describe we set the magnetic field strength to the value B p = 10 15 Gauss. Fig. [T] shows 
the waveform of the gravitational wave observed at r = 2000M, the amplitude is normalized to correspond to an 
emitted energy of £?gw ~ 1-8 x 10 49 (M/50M©) ergs. On the other hand, the waveforms of electromagnetic waves 
driven by the gravitational waves are shown in Fig. [2] From this figure, we can observe somewhat complicated 
waveforms of electromagnetic waves due to the coupling with the gravitational waves. From the specific waveforms, 
one can estimate the value of the proportionality constant in the relation (Eq. (|6.3[0 to be a = 8.02 x 10 -6 . This 
efficiency might not be very high, but the radiated energy of gravitational waves can reach ~ 10 51 ergs for a black 
hole formation due to the merger of a neutron stars binary (see, e.g. [Hj]). In this case the strength of the magnetic 
field can be amplified by the Kelvin- Helmholz instability to reach values of the order of 10 15-17 Gauss (56|. Although 
this is not an ideal situation for the black hole case we are considering in this paper, if one adopts the above efficiency 
for the case of a black hole formed after merger, one can expect that energies of the order of ~ io 46-50 ergs can be 
emitted in the form of electromagnetic waves which can be potentially driven by the gravitational field perturbations. 

Furthermore, in Fig. [31 we show the Fast Fourier Transform (FFT) of the electromagnetic waveforms shown in 
Fig. [5J where for comparison we also add the frequencies of the quasinormal modes for I = 1 electromagnetic waves 
(dashed line) and for I = 2 gravitational waves (dot-dash line) radiated from the Schwarzschild black hole [45| . From 
this figure, one can obviously see that the driven electromagnetic waves have two specific frequencies corresponding 
to the I = 1 quasinormal mode of electromagnetic waves and the I = 2 quasinormal mode of gravitational waves. 
This means that it might be possible to see the effect of gravitational waves via observation of electromagnetic 
waves. However, electromagnetic waves with such a low frequencies could be coupled/absorbed by the interstellar 
medium (and/or accretion disk around the central object) during the propagation and it will be almost impossible to 
directly detect the driven electromagnetic waves. The only possible way to see the driven electromagnetic waves is 
the observation of indirect effects, such as synchrotron radiation. 
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FIG. 1: Gravitational waveform observed at r — 2000M. In the right panel, we also show the absolute value of X20. 




FIG. 2: Waveform of the driven electromagnetic waves observed at r = 2000M for B p = 10 15 Gauss. In the right panel, we 
also show the absolute value of *io. 

B. Perturbations on a Neutron-Star Background 

In order to examine the coupling between the emitted gravitational and electromagnetic waves in a neutron star 
background, we adopt the same initial conditions as for the black hole case, i.e. the electromagnetic perturbations are 
set to zero and the initial gravitational perturbations arc approximated by an ingoing Gaussian wave packet. In the 
numerical calculations, we adopt a grid spacing of Ar = i?/200 and a time step At/ Ar = 0.05, where R is the stellar 
radius. For the background stellar models, we adopt the polytropic equation of state (EOS) of the form P = Kp T . 
Then, one can get the waveforms of the reflected gravitational waves and the induced electromagnetic ones. 




FIG. 3: FFT of the electromagnetic waves shown in Fig. [2] The two vertical lines correspond to the frequencies of quasinormal 
modes for I = 1 electromagnetic waves (dashed line) and for I = 2 gravitational waves (dot-dash line). 
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As an example, we show results for a stellar mode with T = 2 and K = 200 km 2 . Fig. [4] shows the waveforms of 
the gravitational waves (solid line) and the electromagnetic waves (dotted line) observed at r = 300 km, where we 
adopt two stellar models with different compactness M/R (see Table Q] for the stellar properties). Compared with the 
fast damping of gravitational waves, one can see the long-term oscillations in the electromagnetic waves, which can be 
driven not only by the quasinormal ringing of gravitational waves but also during the tail phase of the gravitational 
waves. For the waveforms shown in Fig.21 the FFT is plotted in Fig. where the left and right panels correspond to 
the FFT of the gravitational and electromagnetic waves, respectively. From this figure, one can see the same features 
as in the case of a black hole. Namely, the FFT of the electromagnetic waves driven by the gravitational waves has 
two specific frequencies, i.e. one is the proper electromagnetic oscillation (1st peak in the right panel of Fig. [5]) and 
the other one is the oscillation corresponding to the gravitational waves (2nd peak in the right panel of Fig. [5]). We 
remark that electromagnetic waves with such low frequencies could be absorbed by the interstellar medium and then, 
their direct detection is almost impossible. Namely, we should consider the secondary emission mechanism such as a 
synchrotron radiation. Maybe, the plasma around the central object will be excited after receiving the energy from 
the electromagnetic waves driven by the gravitational waves and move along with the magnetic field lines. Anyway, 
such a secondary emission mechanism will be discussed somewhere. Furthermore, we find that as in the case for a 
black hole, the relationship between the emitted energies of gravitational and electromagnetic waves can be described 
by Eq. (|6.3|) . even for neutron stars, if B p is considered as the magnetic field strength at the stellar pole (r = R and 
6 = 0). In practice, for the specific stellar models in Fig. 21 the proportionality constant becomes a = 1.61 x 10~ 5 
and 4.37 x 10~ 6 for the particular stellar models with M/R = 0.162 and 0.237, respectively. 




t [msec] t [msec] 



FIG. 4: Waveforms of the gravitational waves (solid line) and the electromagnetic waves (dotted line) for the stellar model with 
B p = 10 1 Gauss, which are observed at r = 300 km. The left and right panels are corresponding to different stellar models 
for EOS with F = 2 and K = 200 km 2 . 
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FIG. 5: FFT of the gravitational waves (left panel) and electromagnetic waves (right panel) shown in Fig. [4] 

In order to see the dependence on the stellar properties, we study a variety of stellar models with different stiffness 
of the equation of state and with different central densities, radii, and masses, which are given in Table [I] As a result, 
we find that the proportionality constant a can be written as a function of the stellar compactness, which is almost 
independent of the stellar models and the adopted equation of state. In fact, in Fig. [5] we plot the values of a 
for various stellar models, where the circles, diamonds, and squares correspond to the results for the stellar models 
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TABLE I: Stellar parameters adopted in this article. 
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P< 
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0.194 


2.25 


600 


3 





X 


10 15 


1.404 


8.48 


0.245 


2.25 


600 


4 





X 


10 15 


1.486 


7.95 


0.276 



characterized by (T,K) = (2,100), (2,200), and (2.25,600). From this figure, one can see that the proportionality 
constant a depends strongly on the stellar compactness, as expected, with typical values ranging from 10~ 6 up to 

~ io- 4 . 

lO- 3 
10^ 



B IO" 5 
10" 6 



'"0.10 0.14 0.18 0.22 0.26 0.30 

MIR 

FIG. 6: The proportionality constant a as a function of the stellar compactness for various polytropic models. The circles, 
diamonds, and squares correspond to stellar models with (r, K) = (2, 100), (2, 200), and (2.25, 600). 

VII. CONCLUSION 

We have considered the coupling between gravitational and electromagnetic waves emitted by compact objects, 
i.e. black holes and neutron stars. We have derived a coupled system of equations describing the propagation of 
gravitational and electromagnetic waves. In our study we have investigated the driving of electromagnetic pertur- 
bations via their coupling to the gravitational ones. However, for simplicity, we have neglected the back reaction 
from the electromagnetic waves on the gravitational waves, because the magnetic energy of the compact objects, even 
for magnetars, is quite small as compared with the gravitational energy. We found that the electromagnetic waves 
of specific parity with harmonic indices (l,m) can be coupled to gravitational waves of the same parity and with 
harmonic indices (^, m) (for ra/0) and harmonic indices (I ± 1, m), for every value of m. In particular, our findings 
lead to the result that, for the axisymmetric perturbations, i.e., m = 0, the dipole electric electromagnetic waves will 
be driven by axial quadrupole gravitational waves, while the dipole magnetic electromagnetic waves will be driven by 
polar gravitational waves. 

As an application of our perturbative framework, we presented numerical calculations for the case in which dipolc- 
clcctric electromagnetic waves are driven by the axial gravitational ones, both for the case of a black hole and a 





-•- (2,100) 




-♦- (2,200) 




(2.25,600): 
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neutron star background. Wc found that the emitted energy in electromagnetic waves driven by the gravitational 
waves is proportional to not only the emitted energy in gravitational waves but also to the square of the strength 
of the magnetic field of the central object. For the case of a black hole background, the ratio of the emitted energy 
of the electromagnetic waves to that of the gravitational waves is around 8 x 10~ 6 (-B p /10 15 G) 2 , where B p is the 
magnetic field strength at r = 6M. On the other hand, in the case of a neutron star background, we find that this 
proportionality constant can be written as a function of the stellar compactness. 

Although we have considered only the case of axial gravitational waves and the associated induced electromagnetic 
waves, the polar oscillations also play an important role in extracting the information about the neutron star structure 
since in the case of non-rotating stars, the matter oscillations are typically coupled to the polar gravitational waves. 
This is a direction that we are currently investigating. 
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